arXiv:l509.04968vl [hep-ph] 14 Sep 2015 


Robust characteristics of nongaussian fluctuations from the NJL model 


Jiunn-Wei CherQ 

Department of Physies and Center for Theoretieal Seienees, 

National Taiwan University, Taipei, 10617 Taiwan 
Leung Center for Cosmology and Partiele Astrophysies (LeCosPA), 

National Taiwan University, Taipei, 10617 Taiwan and 
Helmholtz-Institut fur Strahlen- und Kernphysik and Bethe Center for 
Theoretieal Physies, Universitdt Bonn, D-53115 Bonn, Cermany 

Jian Den^ 

Key Laboratory of Partiele Physies and Partiele Irradiation (MCE), 

Sehool of Physies, Shandong University, Jinan 250100, China 

Hiroaki Kohyamsj^ 

Department of Physies and Center for Theoretieal Seienees, 

National Taiwan University, Taipei, 10617 Taiwan 

Lance Labur|i] 

Department of Physies, University of Texas, Austin 
(Dated: 7 September, 2015) 

We evaluate the third- and fourth-order baryon, charge and strangeness susceptibilities near a chi¬ 
ral critical point using the Nambu-Jona-Lasinio model. We identify robust qualitative behaviours of 
the susceptibilities along hypothetical freeze-out lines that agree with previous model studies. Quan¬ 
titatively, baryon number fluctuations are the largest in magnitude and thus offer the strongest signal 
when freeze-out occurs farther away from a critical point. Charge and strangeness susceptibilities 
also diverge at a critical point, but the area where the divergence dominates is smaller, meaning 
freeze-out must occur closer to a critical point for a signal to be visible in these observables. In case 
of strangeness, this is attributable to the relatively large strange quark mass. Plotting the third- 
and fourth-order susceptibilities against each other along the freeze-out line exhibits clearly their 
non-montonicity and robust features. 


I. INTRODUCTION 

Heavy ion collision experiments and lattice simulations 
are probing the phase diagram of QCD matter to help un¬ 
derstand the chiral and deconfinement phase transitions 

[1] . These experiments have shown that at low density, 
the transition is a continuous cross-over at T ^165 MeV 

[ 2 ] . At high density, models lead to the expectation that 
the transition is first order. This structure would be veri¬ 
fied by locating a critical end point in chemical potential- 
temperature plane where the first order line begins. The 
heavy ion collision experiments yield statistical observ¬ 
ables of QCD matter, including proton number and elec¬ 
tric charge fluctuations mil] ,and lattice studies have de¬ 
veloped a variety of approaches to searching for structure 
at /iB 7^ 0 [SUIT]. 

To interpret the experimental data in terms of the 
phase diagram, we need to relate them to signatures de¬ 
rived from theory calculations. Since conventional lat¬ 
tice methods are limited to low baryon density by the 
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sign problem, it is useful to complement the investigation 
with model studies in which we can thoroughly study 
the phase structure at /i 7 ^ 0. We use the models to 
identify robust signatures of the critical point, such as 
consequences of the large correlation length, that would 
be manifest also near a critical point in QCD. The corre¬ 
lation length becomes large near the critical point (CP) 
because the CP is a second-order phase transition where 
the mass term in the Landau effective potential for the 
order parameter a vanishes, rricr ^ ^ 0. Therefore, 

at the CP, the longest wavelength correlations (of order 
the system “size”) can be investigated using the partition 
function [ 12 ] 

2: = (1) 

with the Landau effective action 

Q[cr] = Qo + ^ 

Qo is the energy due to the vacuum expectation value 
(vev) of the order parameter, a is the fluctuation around 
this vev, and for fluctuations at wavelengths similar to 
the system size, we consider only the zero momentum 
mode, so the kinetic term (Vcr)^ vanishes. 

The divergent parts of the fluctuations associated with 
^ ^ oc at the CP (an infrared fixed point) are universal 
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for systems within the same universal class. The finite 
parts of the fluctuations are model dependent which can 
be described by higher dimensional operators in the effec¬ 
tive action. Using Eq. Q, statistical observables such as 
baryon number fluctuations are related to fluctuations of 
the order parameter m- For example, higher order, non- 
Gaussian fluctuation moments are more sensitive to a 
critical point because they diverge with a larger power of 
the correlation length, as shown for the third and fourth 
order susceptibilities [IMS]. 

In this approach, we hypothesize the presence of a 
critical point, investigate the consequences, and compare 
them to data. It is important to bear in mind that the 
conditions we identify for a critical point can only be 
necessary conditions but not sufficient conditions. Also, 
there are many assumptions in our attempt to compare 
with heavy ion data. For example, we assume the fire¬ 
ball is near thermal equilibrium at freeze-out, though 
critical slowing of dynamics would be important if the 
fireball approaches the critical point chemical potential 
and temperature m- Additionally, there can be changes 
in expansion dynamics and interactions that produce 
variations in particle spectra and acceptance indepen¬ 
dent of critical phenomena, whose fluctuations must also 
be controlled for m- To claim evidence of a critical 
point in this context, we should identify as many model- 
independent signatures as possible to corroborate inter¬ 
pretation of the data in terms of a critical point. 

In this work, we study the Nambu-Jona-Lasinio (NJL) 
model, which is a QCD inspired field theoretic model 
with spontaneous breaking of Chiral symmetry HE- 
[20] . The model comprises distinct flavors of quarks 
that interact via (effective) point-like 4-fermion opera¬ 
tors. At low density and temperature, the model ex¬ 
hibits a non-zero chiral condensate {'ip'ip) 7 ^ 0 [2TII25] . 
The goal is to identify, in the NJL model results, char¬ 
acteristics of fluctuation observables that are predicted 
to be model-independent by analysis such as given in 
na ng. Previous work using the NJL model has fo¬ 
cused on third moments m, or electric charge suscep¬ 
tibilities [26]. Numerical studies of the susceptibilities 
have also used the Polyakov-loop improved NJL model 
EH and Dyson-Schwinger approaches [28]. Some of the 
model-independent characteristics have been confirmed 
and cross-checked in the Ising m and Gross-Neveu (GN) 
models [29] . 

There are two main goals in this work. The first goal is 
to make use of the flavor dependence in the NJL model 
to compute the complete set of susceptibilities on the 
phase diagram. The second goal is to check whether 
some mo del-independent features obtained from a gen¬ 
eral effective potential analysis then confirmed by the 
GN model [29] will remain in the NJL model, which be¬ 
longs to the same universality class with QGD at the 
GP. These features include: ( 1 ) The negative cr-kurtosis 
((cr^) — 3(cr^)^ < 0 ) region occurs almost entirely in the 
symmetric (normal) phase. (2) However, in addition to 
the cr-kurtosis, there are more critical-mode correlators 


contributing to the fourth-order baryon number suscepti¬ 
bility \nZ/diJ.%) to determine its negative region. (3) 
The peaks in non-Gaussian susceptibilities on a freeze 
out line obey an ordering in temperature. 

Among those features, we expect ( 1 ) to be robust, be¬ 
cause they can be understood from the shape of the ef¬ 
fective potential around the critical point [29]. We also 
expect ( 2 ) to happen, that is, there are other terms as 
important as the cr-kurtosis to determine the negative 
region of {d^ In Z/d/ j.%). This is based on a general ef¬ 
fective action analysis of Ref. [29] to identify terms with 
leading power divergence in ^ near the GP. We will check 
it numerically using the NJL model. 


II. COMPLETE SET OF FLUCTUATION 
OBSERVABLES 

Fluctuations in conserved charges are important ob¬ 
servables because they can be obtained in principle from 
HIC experiments as well as lattice simulations. Not all 
fluctuation moments can be measured in practice, and 
we will discuss below the observationally most relevant 
subset. The fluctuation moments are derivatives of the 
partition function with respect to the chemical potentials 
of the conserved charges: 

5 5 ••• \^) 

li a = f3 = .... we shall also use = d^lnZ/dfi ^. 
Recall the first derivative is the expected number, 
d\nZjdjix = {Nx)^ which is conserved and there¬ 
fore constrained by the initial conditions. For example, 
since protons and neutrons comprise only up and down 
quarks and strangeness is conserved by QGD reactions, 
{Ns) = 0. On the other hand, since the heavy nuclei 
typically contain more neutrons than protons, the initial 
state has an isospin asymmetry. 

At quark level, we have three independent chemical 
potentials jiu^ jidi and jig associated with the conserved 
quark numbers for d, and s quarks, respectively. To 
compare to experiment, we use the basis of conserved 
charges {B^Q^S}. The hadron-level strangeness (large 
S) is defined so that has S' = -hi in agreement 
with the experimental notation and opposite the quark- 
level definition of the strange charge (small s) where the 
strange quark has s = —1 and antistrange quark has 
5 = +1. We will use the {g, /, s} basis of charges, where 
2iiq = fjiu ^ lid and iij = fiu — fJ^d- Since the isospin po¬ 
tential III is also measurable at hadron-level, we can use 
observations to constrain the chemical potentials, and so 
also the set of susceptibilities we need to evaluate. 

The chemical potentials MS',at chemical 

freeze-out are determined by fitting the observed par¬ 
ticles yields to a statistical model of hadronization in 
the fireball. Limited acceptance and experimental effects 
can mean the fit values of the chemical potentials differ 
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from the theoretical expectation. The data indicate that 
when collision energy per nucleon pair in the center-of- 
mass frame {^/sJ^) is 7.7, 39, 200 GeV, is ^130, 
30, 10 GeV and jj^s is ^90, 20, 5, respectively [30]. Using 
Eq. ^ and isospin symmetry /i^, this implies fig is 
^40, 10 , 5 GeV, even though theoretically {Ng) = 0 in 
the initial state. However, given jiig <40 MeV which is 
smaller than the lower range on the strange quark (cur¬ 
rent) mass rrig > 90 MeV, it is a good approximation to 
set fig = 0 in our model calculations: dependence of the 
phase diagram on fig is weak as long as for fig < rrig. 

It is worth noting here that one piece of evidence 
for non-equilibrium at chemical freeze-out is that par¬ 
ticle yield fits are improved by including a strange quark 
quenching factor 7 ^, which observationally is < 1. This 
smaller-than-equilibrium abundance of both s and s can 
be understood considering that the phase space for pro¬ 
ducing s, s quarks is smaller due to their moderately 
large mass rrig > 90 MeV compared to the temperature 
T < 160 MeV. That is, a substantially smaller fraction 
of gluons (or qq pairs) collide with center-of-mass energy 
> 2mg > 180 MeV, which suppresses the reaction rate 
for production in comparison to the lifetime of the fire¬ 
ball [31]. Although it impacts the absolute abundance 
of measured strange particles, it is linearly independent 
of fig that controls net strangeness, i.e. the number of 
strange minus the number of antistrange quarks Ng—Ng. 
With /ig = 0, there is no violation of strange quark num¬ 
ber independent the value of 7 ^. 

Additionally, analysis of the charged pion ratio 
V( 7 r“)/V( 7 r<) has consistently concluded that fij = 0 
at high beam energies, though fij can differ from zero at 
the lowest collision energy ^/S]sfN = 10 GeV [32]. This 
variation as a function of can be understood by 

the total particle multiplicity becoming much higher for 
larger ^/snn- We have checked quantitatively in the NJL 
model, that the phase diagram is weakly dependent on 
fij as long as fij <C 65 MeV (this is also found in 

[33]b Therefore, in our work we shall for simplicity set 
fii = 0 for all fiB- 

With /ig = yU/ = 0, a complete set of nontrivial suscep¬ 
tibilities up to fourth order is: 

XqqiXlhXss^ 

xf\Xqss,XqII, 

, xf^ ,XlIss, 

These are easy to count in the {g, /, s} basis, because 
/^/ = = 0 implies that odd-order derivatives of fij 

vanish identically. When we choose another basis, there 
can only be 3, 3, and 6 independent susceptibilities at 
2nd, 3rd and 4th order respectively. These are plotted in 
Figs. [^ [^ [^ and 

While the long-wavelength correlations are dominated 
by the iso-scalar critical mode a and the isospin chem¬ 
ical potential is “small” (as defined below) this relation 
shows that charge fluctuations diverge like the baryon 
fluctuations near the critical point [34]. As explained 


in the next section and Appendix the second-order 
isospin and strange susceptibilities do not have a singu¬ 
larity at the critical point, even in the presence of flavor 
mixing. By direct computation of XQ Xb in the 
NJL model below, we will verify this fact. On the other 
hand, because the correlation length in the fireball may 
peak at ^ 3 fm, only about twice the thermal correlation 
length (and much smaller than infinity), the difference 
between the singular iso-scalar critical mode contribu¬ 
tion and the smooth model-dependent contributions to 
correlations may be less dramatic. 

To relate the B^Q^S susceptibilities to the derivatives 
with respect to quark susceptibilities we need to rewrite 
partition function in terms of the B^Q^S charges. The 
grand canonical partition function can be expressed as 

2 : = (4) 

where Nx is the number operator, i.e. the time com¬ 
ponent of the conserved current. This implies that the 
quark chemical potentials fiq^fij^fig are related to the 
chemical potentials associated with observable numbers 
fiB^fJiQ^fis using the linear relations between the particle 
numbers. Starting from the definitions 

Ng = N^ + Na, Ni = ^{N^-Na), 

we have 

NB = l{Nq + N,) = ^{Nu + Nd + N,) (5) 

NQ = lNg + Ni-^N, = ^ (27V„ -Na-N,) (6) 

Ns = -N, (7) 

These relations lead to 

IJ'B = 3/ig - (8a) 

Hq = Hi (8b) 

M5 = - Ms (8c) 

and conversely 

11 

Mg = + -MQ (9a) 

M/ = MQ (9b) 

11 

Ms = gMs - gMQ - Ms (9c) 

Using Eq. § we can rewrite the observable susceptibil¬ 
ities in terms of quark susceptibilities, which are deriva¬ 
tives of the NJL potential. For example, the baryon sus¬ 
ceptibility is 

d^lnZ (Id 1 5 ^ 

5m| “ VsS/xUsa/xsj "" 

12 1 

— gXqq T ^ 







Since net strangeness is conserved and equal to zero 
(Ns) = 0 independent of the d chemical potentials, 
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( 11 ) 


SO that the second term in Eq. (10) vanishes. Finally, we 
have 


Xb 

Similarly, the charge susceptibility simplifies 

XQ = (13) 

because the cross term (x Xus is zero. Note the pres¬ 
ence of the strange susceptibility. This shows that both 
the charge and baryon fluctuations are dominated by the 
light quark fluctuations, as we expect [34]. The strange- 
hadron susceptibility is 


= Xs 


(14) 


III. SUSCEPTIBILITIES WITH FLAVOR 

In Ref. [29|, the explicit expressions for the quark 

(baryon) susceptibilities Xb^ ? Xb^ terms of autocor¬ 
relation functions of the a zero mode {a'^) were derived 
for the GN model. The k-th order susceptibility requires 
only correlation functions of k-th order and lower, see 
Eqs. (6), (7), and (8) in [29|. While higher order sus¬ 
ceptibilities typically diverge with larger powers of the 
correlation length ^ m, we found that the correspond¬ 
ing autocorrelation function {a^) does not always pro¬ 
vide the complete leading order contribution to Xn- H 
also implies that the evaluation of Xn can be organized 
diagrammatically corresponding to a finite set of terms in 
perturbation theory. As we extend the analysis to NJL, 
the diagrammatic expansion becomes more complicated, 
however, we can use it to help explain the behaviour of 
strange quark susceptibilities near the critical point. 

When writing these susceptibilities in terms of order 
parameter fluctuations, we must include three conden¬ 
sates, one for each flavor, and hence three fluctuation 
fields ctq,, (a = u^d^s. Due to the anomaly-induced six- 
quark interaction vertex, the a “correlator” (statistical 
covariance) &^VtjdGoidG^ is not diagonal in flavor. This 
fact is derived and a brief description of the diagrammatic 
system is found in Appendix [Bj 

We have checked numerically that the expansion in 
terms of g correlators is equivalent to taking the /i deriva¬ 
tives directly at the minimum of the effective potential 
by brute force, e.g. a finite difference method. The g- 
correlator diagrams provide an easier method to organize 
and analyze the many terms arising from the evaluation 
of higher order and multi-flavor susceptibilities. 


O 

FIG. 1. Diagrams (a) and (b) resp ectively correspond to the 
first and second terms in Eq. ( |B3| . Each crossed circle de¬ 
notes an insertion of the q^qa operator arising from taking a 
fia derivative to the pressure which is proportional to the log¬ 
arithm of the partition function. The solid and dashed lines 
denote the quark and the a a propagators. The cross denotes 
the flavor mixing of the a propagators. 


Using the diagrammatic system developed in Appendix 
[B] we can quickly identify large contributions to flavor- 
dependent observables. As mentioned above, the critical 
mode is associated with light quarks, since this transi¬ 
tion occurs at smaller /r^. Therefore the most divergent 
contributions are n-point correlation functions of the Gq 
fluctuations. 

The diagrams shown in Figure [^explain the behaviour 
of the second order susceptibilities xii Since the 

first derivatives with respect to fij and fis vanish iden¬ 
tically at = 0, the diagram on the right vanishes 

because it has only a single /i-derivative acting on the 
distribution function. As a result, there are no contribu¬ 
tions involving cr-correlators to these susceptibilities, and 
they show no divergence. 

In Figures iilD and we display the complete set 
of quark susceptibilities Xa/3... up to the fourth deriva¬ 
tives, from which we need to construct the fluctuations 
of observables. 

The third order susceptibilities are straightforwardly 
understood since each is a single fiq derivative on a second 
order susceptibilities. Thus, we can confirm by eye the 
correct behaviours: Xqqq is odd across the first order line, 
corresponding to the peak seen in Xqq- Xqii hus a peak 
around the first-order line, corresponding to the jump in 
Xii but remains positive at high /i since xii continues 
with positive derivative. The variation of Xss, exhibited 
by nonzero Xqss^ is too small in magnitude to see in Fig. 
1^ The derivative with fiq contributes a singularity at 
the GEP to each Xqii und Xqss, which arises from a new 
nonzero diagram, shown in Fig. 

Among the fourth order susceptibilities, three 
X^\xqqii^Xqqss cau be uuderstood in the same way 
as the third order susceptibilities looking at the 
derNatives of the previous plots. For the other three, 
we can use the diagrammatic analysis. 

Xs^^ can have a singular contribution. Again, the dia¬ 
grams contributing are severely constrained by the fact 
that odd derivatives in fig vanish. The only two dia¬ 
grams with an even of number of fig derivatives on the 
external legs are shown in Figure The diagram on 
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FIG. 2. Xqq^Xii^Xss, three independent components of the 
second order susceptibilities for the conserved charge {q, /, s). 


the right does contain a a correlator, which is flavor- 
diagonal, connecting one 5 distribution function (loop) 
to a second. This correlator is also singular at the light 
quark critical point, due to the flavor coupling term in 
the NJL Lagrangian. This fact is recognized in many 
studies of NJL m, that the strange condensate also has 
a small discontinuity across the first order line. However, 
the singularity at the CEP is suppressed quantitatively 
by the small discontinuity and large strange quark bare 
mass. With high precision numerics, it becomes visible 
in Figure pT| below. 

The same diagrams in Fig. H contribute to With 
the presence of cr-correlators for the light quarks, we can 





FIG. 3. Txqqq,Xqii^Xqss, thrcc independent components 
of the third order susceptibilities for the conserved charge 




FIG. 4. Diagram contributing the singularity to Xqii 

Xqss- 
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FIG. 5. The fourth-order susceptibilities, ^ind Xs^^- 

expect some singularity, larger than in Xs^^ but smaller 
than in x^^ bue to some cancellation of leading terms. 
For Xiiss only the diagram on the right in Fig. con¬ 
tributes (because the {qiqi)^ {qsqs)‘^ operators cannot be 
contracted to a single loop). The cr-correlator implies 
Xiiss is singular, but again the divergence is relevant only 
in a very small region, even accounting for the ^-scaling of 
the coefficients weighting each a-correlator contribution. 

The cross-correlation in flavor Xqqss receives a larger 
singular contribution. To see this, consider the four pos¬ 
sible diagrams shown in Figure In each diagram, there 
are at least two (jq propagators, which together give a 



FIG. 6. The fourth-order susceptibilities, Xgqii^Xqqss and 

Xiiss- 
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O D -0 

FIG. 7. Diagrams in the expansion of Due to the /Xs 

—fis symmetry, the only nonvanishing diagram involving a a 
field correlator is at right, and can only involve the as order 
parameter field. 
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4- 


'Q 


FIG. 8. Diagrams in the expansion of Xuuss 


factor oc Further, we can estimate that the dominant 
effect near the CEP is the diagram with the au 3-point 
function: involving 3 a correlators (each scaling as 
and the 3-point vertex (scaling as (^T)“^-^, obtained by 
solving the gap equation near the CP), the diagrams has 
an overall scaling ^ before the ^-dependence of the 
coefficient is applied. The dominance of this term is con¬ 
firmed numerically. However, the overall magnitude of 
the peak is exponentially suppressed by the strange mass 
as compared to the singularities in light quark suscepti¬ 
bilities, which we have checked by varying the bare mass 
of strange quark. 


Writing out the nongaussian hadron-level susceptibil¬ 
ities, as expected, only 3 third-order susceptibilities and 
6 fourth-order susceptibilities are linearly independent. 


The third order susceptibilities are 


(3) 

Xb 


JLy(3) I ly 
(2IJ Xq ' 


= d_y(3) 
216^^ 

xf = 0 


1 


-)XqII 


Xbbq = 
Xbbs = 
Xbqq = 


_ —y 

Xq Xqs 


54" 


:Xqs 


1 


108 


18 


Xqs 


Xqqq ' ^Xqll 

1 


1 

^XBSS — XBQS — ^XQSS — -jXQSS — Y^Xqs 


and the fourth order susceptibilities are 


(15a) 

(15b) 

(15c) 

(15d) 

(15e) 

(15f) 

(15g) 
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2 
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+ ^^Xqq. 
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2 
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1 
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- Xy(4) 

27AS 

Xbbqq 

- d_y(4) . 
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Xbbss 

Xbqqq 
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Xqqqs 
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+ gX//SS + ^xl^^ 

■ ly + 

^Kqqss 4^ gYs 
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1 


gX/Zss - ^xl^^ 
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_ _i _1 (4) 

1 _1 (4) 

-^^Xqqss gXs 

-^Xqqss + Xliss + ^Xs 

1 1 74 ) 

'^Xqqss ffi Xliss T 

1 ] 


1 1 (4) 

-^Xbsss = -^Xqsss = 


(16a) 

(16b) 

(16c) 

(16d) 

(16e) 

(16f) 

(16g) 

(16h) 

(16i) 

(16j) 

(16k) 

(161) 

(16m) 


Based on the above analysis, the behaviour of these sus¬ 
ceptibilities near a critical point (magnitude and shape) 
are dominated by the light quark susceptibilities, such 
as x^q^ third order and x^q^ fourth order. Only in 
observables without these two terms can the singularity 
due to the strange quark fluctuations be visible. We will 
see these effects in the following section. 




FIG. 9. mi for S, Q on the phase diagram, mi for S' is zero. 


IV. CHARACTERISTICS OF THE 
SUSCEPTIBILITIES 


Here we display the numerical results for the fluctua¬ 
tion observables of primary interest. The volume factors 
in the susceptibilities are removed by considering ratios 


mi{X) 


Txf ^ 

Xxx ’ 


m2{X) 




(17) 



Note that 


explained 


for X = B^Q^S. We display the results for mi,m 2 for 
each conserved charge in Figures and 
mi (S') oc ^555 = 0, so we do not show it. 

In m 2 (S), we see the divergence in 
above. However, as seen in the scaling of the axes, 
the magnitude of the singularity and area of the crit¬ 
ical region are small. It requires high numerical pre¬ 
cision to make the effect visible. Similarly we can 
form observable ratios from other third- and fourth- 
order hadron-level susceptibilities, involving strangeness 
or cross-correlations in baryon and charge fluctuations. 
The ratio Txbss/xbb is proportional to Xqss: however, 
the divergence in Xqss is effectively canceled by the diver¬ 
gence in Xqq^ ^s suggested by comparing Figand Figj^ 
Consequently, the ratio shows no singular behaviour at 



FIG. 10. m 2 for each of S on the phase diagram. 
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the CEP. For the ratio T'^Xbbss/xbb, Xuuss in the nu- 
merator diverges with a larger power of ^ than the de¬ 
nominator, as seen by comparing the diagrams in Figj^ 
and Fig{^ so there is a visible singularity peak at CEP. 
The shape of T‘^Xbbqq/xbb is very similar to m 2 (5) 
and m 2 (Q), because all three are dominated by the light 
quark susceptibilities. The results in Fig. [TT] confirm that 
singularities in the isospin and strangeness susceptibili¬ 
ties make small impact on the observable fluctuations. 

The area of the positive and negative regions in the 
mi, m 2 for 5, Q are seen more clearly in the density plots 
Figs, an d The shapes are consistent with previ¬ 
ous results |14j. Note the area of the critical region for 
the baryon number fluctuations is larger than the same 
region for charge fluctuations, consistent with the overall 
larger magnitude of the baryon fluctuations seen in Figs. 

and The area of the critical region for strangeness 
is smaller still than for the charge fluctuations, and is 
due to the small effect of the singularity on the strange 
quark susceptibilities away from the critical point. In 
the context of the conventional equilibrium fluctuations 
framework to look for signatures of a critical point, the 
large mass of the strange quark (comparable to T,/ig) 
suppresses fluctuation observables. However, this analy¬ 
sis does not account for the possible impact of nonequi- 
librium dynamics, such as critical slowing, on strangeness 
yields and fluctuations. 

In addition, following na and [29] , we show the result 
of extracting the values of mi and m 2 along various hypo¬ 
thetical freeze-out lines chosen to pass varying distances 
from the critical point. The freeze-out lines are shown in 
the top frame of Fig. 14 The behaviour of mi{B) and 


m 2 {B) along the freeze-out lines shows good qualitative 
agreement with the previous GN model calculations, sup¬ 
porting the robustness of these shapes and their variation 
with distance from the CEP. The profile of charge fluc¬ 
tuations along the freeze-out line is very similar to the 
baryon fluctuations, with the apparent difference due to 
the overall smaller amplitude of the variation. In Eigflbl 
we show the ratios involving strangeness along the freeze- 
out lines, Bxbss/xbb, ^ 2 ( 6 ') and Bxbbss/xbb- As 
the critical region for strangeness is very small, these ra¬ 
tios are not sensitive to the singularity of the strange 
quark susceptibilities. 

Einally, we show how mi and m 2 vary in relation to 
each other along the freeze-out lineQ In Eig. [T^ de¬ 
creasing x/snn (increasing fiB or decreasing T) starts 
from mi = 0 and continues in an anti-clockwise trajec¬ 
tory around the loop. The convergence of the lines at 
high T (close to mi = 0) is due to the convergence of the 
hypothetical freeze-out lines in our modeling, see Eig. 
(top). Eor the lower T freeze-out lines (green dashed 
and blue dotted lines), the loop does not close. Even so, 
the low T /high ji end of the trajectory is the same for 



. .■■ 

T (MeV) 




T (MeV) 
45 // 



FIG. 11. Ratios of cross-correlations T^^XbbqqIxbb, 
T‘^Xbss/xbb, and T‘^Xbbss/xbb on the phase diagram. 


^ We thank R. Gavai for suggesting this plot. 


















10 


200 


150 


QJ 


100 


50 


200 


150 


QJ 


100 


50 


mi(B) 



100 200 300 400 

Hq (MeV) 


m2(B) 



100 200 300 400 

fiq (MeV) 


200 


150 


QJ 


100 


50 



i 


200 


150 


100 


50 




I 


0 100 200 300 400 

Hq (MeV) 


miiQ) 



0 100 200 300 400 


(MeV) 


FIG. 12. For baryon susceptibility, the figures show the sign 
change of mi and m 2 on the phase diagram. Red areas are 
negative, and blue are positive. The green point indicates the 
location of CEP. The dashed thin line and the solid thick line 
stand for the crossover line and the hrst order phase transition 
line respectively. 


all freeze-out lines because the statistics and fluctuations 
are given by m, /i-dominated thermal distributions there. 
These plots show clearly the ordering of features 

^min,m2 ^ ^max,mi ^ ^max,m2 ^ ^CEP 

which was found also in the GN and Ising models. (In 
the Ising model, we plot the kurtosis of the magnetization 
d^M/dH^ versus the skewness d^MjdH^ along lines of 
constant yielding a plot of the same shape as Fig. [T7| ) 
Again, the magnitudes of the mi, m 2 are much greater for 
baryon fluctuations than for charge fluctuations. We ex¬ 
pect these characteristics are robust. Nonequilibrium ef¬ 
fects and higher-order, model-dependent corrections will 
perturb the loop. However, given also the similarity to 
the GN and Ising results, it appears large corrections 


FIG. 13. Similar to the previous hgure, but for electrical 
charge susceptibility. 


would be necessary to affect the ordering in Eq. (18). For 
ratios with strangeness, this ordering is not seen, because 
the critical region for strangeness is too small to impact 
the observables along the chosen freeze-out lines. 


V. CONCLUSION 

In this paper, we have investigated the behaviour of 
non-Gaussian fluctuations of baryon number, electric 
charge and strangeness on the phase diagram. Our pur¬ 
pose was to characterize their dependence at /i 7 ^ 0 and 
find robust qualitative features that can be used to help 
interpret experimental data in the search for a QCD crit¬ 
ical point. To this end, we evaluated the third- and 
fourth-order susceptibilities in the NJL model, which is in 
the same universality class as QCD. Our analysis is lim¬ 
ited to the mean-field approximation which means that 
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T (MeV) 



T (MeV) 


FIG. 14. Top frame: Hypothetical chemical freeze-out lines on 
the phase diagram to model the behaviour of the fluctuation 
ratios as a function of • Middle and bottom frames: mi 

and 7712 of baryon number susceptibilities along these freeze- 
out lines, identified with lines in the top frame by color and 
dashing. 


we have not included several known corrections: modi¬ 
fications to the critical exponents controlling the diver¬ 
gence of the correlation length at the CEP, and model- 
dependent corrections described by higher-dimensional 
operators in the low energy effective field theory. The cor¬ 
rections to the critical exponents are better known quan¬ 
titatively, but model-dependent corrections are less well- 
controlled and their importance is not known in the con¬ 
text of the physical fireball where the correlation length 
of the critical mode may not be much larger than thermal 
correlation length. 

We have found several characteristics of the third- and 



T (MeV) 



T (MeV) 


FIG. 15. Electrical charge mi and m 2 along the freeze-out 
lines defined in the top frame of Figure 


fourth-order susceptibilities that are likely to be robust. 
These are: 


1. Baryon number fluctuations are the largest in am¬ 
plitude. This is because they are dominantly com¬ 
posed of the leading singular contributions from 
the light quark susceptibilities. Charge fluctua¬ 
tions are suppressed by numerical coefflcients, as 
seen by comparing the expansions in Eq. (15) and 
Eq. (16). Strangeness fluctuations are the smallest 
in amplitude, though at high order they also show 
a singularity at the critical point. The strangeness 
fluctuations are suppressed by the relatively large 
bare mass of the strange quark, rris ^ T. 


2. The “critical region” where we see non¬ 
monotonicity in mi and m 2 along hypothetical 
freeze-out linesis largest for baryon number fluc¬ 
tuations. Therefore it provides the largest signal 
for freeze-outs occurring farther in the /i, T plane 
from the CP. This act follows from the amplitude 
of the fluctuations being largest. 

3. The qualitative profile of the third and fourth order 
susceptibilities along the freeze-out line is model- 
independent. mi for both 5, Q has a single peak 
at a temperature greater than the CP. m 2 has first 
a minimum and then a peak as the critical point is 
approached from high T (equivalently larger colli¬ 
sion energy -\/snn)- The m 2 < 0 negative region 
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T (MeV) 
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FIG. 16. The strangeness-dependent ratios Txbss/xbb, 
7712 ( 6 ') and T‘^Xbbss/xbb along the chemical freeze-out lines 
defined in the top frame of Figure p!4| 


may be accessible at chemical freeze-out in the sym¬ 
metry broken phase. These behaviours have been 
seen in the GN and NJL models. 



T Xbss/Xbb 

FIG. 17. Top and middle: m 2 plotted against mi along the 
chemical freeze-out lines, for both baryon number and elec¬ 
tric charge susceptibilities. Bottom: an analogous ratio for 
the cross-correlation between baryon number and strangeness. 
Golor and dashing of the line identifies it with the freeze-out 
lines defined in Fig.p^ 


Appendix A: NJL model 


4. The features of the rui and m 2 freeze-out profiles 
obey a numerical ordering temperature or collision 
energy, given by Eq. (18). This is demonstrated by 
plotting 7772 versus mi along the freeze-out line, and 
we have argued this trajectory is difficult to change 
in qualitative features. 


It is our hope that these features can aid the interpreta¬ 
tion of experimental data for the third- and fourth-order 
susceptibilities. Although our results may support pre¬ 
liminary indications of critical behaviour in the data, it 
is clear that much more work is required to identify these 
signatures with a QCD critical point. 


In our study, we consider the 3-flavor NJL model with 
two degenerate light quarks and a heavier strange quark. 
The quark level Lagrangian is 

7= T] + GTr[{tpXatp)^ + 

i=u,d,s 

- K(det 5 ( 1 - 75)9 +detg(l+ 75 ) 9 ) (Al) 

The model parameters are the “bare” quark masses mi 
and the 4-fermion couplings (7, K. In addition, for nu¬ 
merical evaluation of the effective potential below we 
must introduce a momentum cutoff A. These model pa¬ 
rameters are fixed by matching the pion, kaon, 77 ' masses 
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and the pion decay constant in vacuum, = 138 MeV, 
rriK = 495 MeV, m^/ = 958 MeV and = 93 MeV. We 
use the residual freedom to set the light quark mass near 
the physical point, so that the parameter set is 

rriu = rrid = 5 MeV, rris = 136 MeV, (A2) 
A = 631 MeV, GK^ = 1.83, KK^ = 9.29 (A3) 


The phase diagram is determined by solving the cou¬ 
pled set of gap equations 



2K E 


(A8) 


Spontaneous symmetry breaking leading to nonzero light 
quark condensates {uu ), {dd) ^ 0 has greater impact 
on phenomenology than the explicit breaking from their 
masses, as is the case in QCD. 

To solve the model and study the phase diagram, we 
take the mean field approximation. After shifting the 
quark bilinears with their vacuum expectation values 
qq {qq) + (pq, the quark fields can be integrated out, 
and we obtain an effective potential depending only on 
T, fiq and ipq {q = d, s). As long as the isospin chemical 
potential 

777 

^65 MeV (A4) 

the 14, d condensates are equal and transition at the same 
/i,T. Pseudoscalar, vector and pseudovector diquark 
condensates are also suppressed and can be ignored. Ex¬ 
perimental data indicates that this condition holds even 
at the lowest collision energies. The resulting effective 
potential has the form 


taking the solution that corresponds to the global min¬ 
imum of the effective potential, (lea^yol tensor is +1 in 
entries with a ^ ^ ^ 0.) We solve the system using two 
independent numerical methods, as a quantitative check 
on our results. 


Appendix B: Diagrammatic system 

/i-derivatives of the pressure include a term for each 
flavor 

dPjT, jl, a) ^ ap ap 

diia dua ^ dap -( 0 ) dfia 

The vectors /I, a are shorthand for 
cTg), and subscript indicates the 
derivative is evaluated at the minimum of the potential. 
The last factor dcrp/dpa can be rewritten using the fact 
that the gap equation is independent of the chemical 
potential 


Q, = 2G{al + al + a^) - AKauadas + y] 12/ (A5) 

f=u,d,s 

% = -2N, I ^ (S/0(A2 - y) (A6) 

+ Tln(l + + Tln(l + 

where Q{x) is the Heaviside step function, Ej = + 

The last two terms in Qf are respectively the particle 
and antiparticle fermi distributions including the chemi¬ 
cal potential. The Mf are effective masses, functions of 
the condensates 


Mf=mf- AGcjf + 2Kcjpcjf., {f ^ f'^ f") (A7) 

The term off-diagonal in flavor implies that the strange 
quark condensate is discontinuous where the light quark 
condensates are discontinuous and vice versa. The larger 
bare mass of the strange quark means the impact of ex¬ 
plicit chiral symmetry breaking is larger than for the light 
quarks, and its phase transition occurs only for larger 
chemical potential and temperature, see [29] for discus¬ 
sion how the position of the critical point depends on the 
bare quark mass. If a critical point is accessible to lattice 
or HIC, it will be the light quark critical line, which is at 
lower chemical potential and temperature. However, due 
to the flavor coupling, some evidence of the criticality 
will be manifest in the strange susceptibilities. We will 
discuss this point further below. 


0 = 


d / dP 
dua \dap 
d^P 


r(o) 


(B2) 


duadap 


f ( 0 ) 


E 


d‘^p 


dapda^ 


r(0) 


dpa 


Note that the second term contains a factor which has 
the form of a two-point correlator of the a field. 

When we write out a general (second order) suscepti¬ 
bility, we find two terms 


d‘^P 

dpidpj 

d‘^p 


E 


d^P 

dfiada-j 


d^P \ ^ d^P 
da^das gm) dapd/ip 
(B3) 


A similar relation has been discussed by [35]. The first 
is the second derivative at the minimum of the potential, 
and the second relates /i-derivatives at different points 
through the a two-point function. We express this natu¬ 
rally by the diagrams in Figure This “a correlator” 


da^dap 



2Ar 


(^e 


dA^ 

dap 

(B4) 


OA^ _ _ 

dap drni^ 


4G6p^ - 2K E 
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is not diagonal in flavor space due to the anomaly- 
induced interaction (proportional to K). 

The third and fourth order susceptibilities include new 
terms, such as the a three- and four-point functions. 
By writing out all the diagrams constructed from these 
pieces, one may search higher order susceptibilities for 
their most singular contributions. 
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